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ABSTRACT 

Experimental data on the stress-strain behavior of a polymer-multiwall 
carbon nanotube (MWCNT) yam composite are used to motivate an initial study in 
multi-scale modeling of strength and stiffness. Atomistic and continuum length 
scale modeling methods are outlined to illustrate the range of parameters required 
to accurately model behavior. The carbon nanotubes yams are four-ply, twisted, 
and combined with an elastomer to form a single-layer, unidirectional composite. 
Due to this textile structure, the yam is a complicated system of unique geometric 
relationships subjected to combined loads. Experimental data illustrate the local 
failure modes induced by static, tensile tests. Key structure-property relationships 
are highlighted at each length scale indicating opportunities for parametric studies 
to assist the selection of advantageous material development and manufacturing 
methods. 

INTRODUCTION 

Since their discovery in 1991 [2] carbon nanotubes have sparked the 
imagination of many research groups around the world [3], Carbon nanotubes can 
be found in several distinct formations including single-wall, double-wall, and 
multi-wall [4], Recently, Baughman’s group at University of Texas has developed a 
method for spinning carbon multi-wall nanotubes (MWCNTs) into a continuous 
yam stmcture. This process is amenable to standard textile manufacturing 
techniques and allows for variation in the number of strands, number of plys, and 
geometric factors such as twist and yam diameter [5], [6], It is the objective of this 
paper to consider a yam stmcture and to explore the atomistic and continuum 
relationships available to describe the elastic modulus and initial strength of this 
novel material system. This paper will highlight several key modeling challenges 
and outline the parameters of interest for textile modeling. Experimental data from 
recent tests by the authors will be presented to illustrate the similarity of MWNT 
yams and natural fiber yams. 
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As presented in Gates [7], traditionally, research institutions have relied on 
a discipline-oriented approach to development and design with new materials. It is 
recognized, however, that within the scope of materials and structures research, the 
breadth of length scales may range more than twelve orders of magnitude, and 
different scientific and engineering disciplines are involved at each level. 
Investigation of the mechanical response of the MWCNT yam material may require 
material level descriptions that span at least six orders of length scale. As outlined 
in this paper, the yam is a stmctural construct that exhibits properties that are 
functions of precursor material, loading conditions, and manufacturing induced 
helical geometry. 

Most approaches to multiscale modeling are viewed as either bottoms up or 
top-down. Ultimately, the objectives of these two approaches, providing predictive 
methods for bulk behavior and guiding materials synthesis, are the same, but the 
methods, techniques, and the cross-scale handshaking is dictated by the approach. 
Ideally, both approaches should provide a two-way path for passing intrinsic, 
critical information across the scales. It is the former approach, bottoms-up, that 
will be the focus of this paper. The underlying assumption is that the intrinsic, 
critical structure and constitutive relationships developed at each scale are passed to 
the next scale with a well defined handshaking method. One of the difficulties 
associated with this assumption is that the information passed through the 
handshake should retain the detail and complexity needed to ensure accuracy. The 
approach considered in this paper requires a handshake but does not require direct 
coupling of computational methods. This approach contrasts with explicitly 
coupled models where full atomistic or molecular scale information is retained in a 
few, preassigned critical regions of the material, and higher scale continuum 
models are then used to model deformation in regions of the material that are more 
remote from the atomistic region [8], 

MATERIALS 

The subject material test at NASA Langley Research Center is a yam- 
elastomer composite. The composite is composed of a silicone rubber (Silicones, 
Inc GI1000) polymer with a single layer of embedded, aligned, four-ply MWCNT 
yam fibers. Figure 1 provides a photograph of a cross-section of a typical test 
coupon. 




Figure 1. a)Optical micrograph of a typical uniaxial tensile test coupon gage cross-section, b) 
Inset of right-most yarn on edge, c) Pictorial representation of area fraction of inset (~75% 
polymer). 

The precursor to the yam structure is the MWCNT. The yam material was 
supplied by Professor Ray Baughman, University of Texas at Dallas. The MWCNT 
yam is created by simultaneously drawing and clockwise twisting (20 tums/mm) 
vertically aligned MWCNTs from a MWCNT forest to form a single yarn. Four of 
these singles are then twisted counterclockwise (4 tums/mm) forming a 4-ply yam 
(Figure 2). The MWCNT yam is similar in many respects to standard textile 
materials [9], and can be viewed as a unique structure whose properties are 
dependent not only on the precursor material (MWCNTs) but on yarn geometry 
and manufacturing methods. Additional details on the yam and the associated 
properties can be found in [10], [1 1], and [12], 
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Figure 2. SEM micrograph of a longitudinal view of typical four-ply multi-wall carbon 
nanotube yarn (2500x; 15kV). 




Figure 3. SEM photomicrograph of a hybrid glass fiber 3-D braid with CNT yarns visible as dark 
lines (left picture), and a 3-D woven fabric with 25-ply CNT yarns incorporated as the first five Z- 
yarns from the bottom (right figure). Ref. Bradford [1]. 

The MWCNT yam can also be used in higher complexity forms such as the 
3-D braid and woven fabric shown in Figure 3. In terms of length scales for the 
yam structure under consideration, the following approximate structure-length 
relationships are assumed: MWCNT(10~ 9 m), fibers(10~ 7 m), single yam(10" 6 m), 
yam plies(10' 3 m). In this section of the paper, both atomistic and continuum 
mechanics approaches will be considered. The intent is to provide an initial review 
of methods available to calculate critical parameters for prediction of stiffness and 
strength. 

Molecular Level Simulations 

The basic unit of the yam is the MWCNT. In order to investigate the 
deformation and forces on a bundle of nanotubes, molecular dynamics simulations 
were performed of a system of (10,10) single-walled carbon nanotubes. Each 
nanotube contained 880 carbon atoms per 5.425 nm length. The nanotubes were 
simulated in a periodic box which contained 9 nanotubes (Figure 4). Each nanotube 
was represented with the Brenner-Tersoff potential [13]. The intra-nanotube 
interactions were modeled with the Lennard-Jones potential. The box was subjected 
to compression in the transverse x and y- directions. No change was made to the z- 
dimensions. The percent of deformation relative to the starting configuration is 
included in Table 1. Altogether, five levels of compression were simulated up to 
15%. 

At each level of compression strain, two simulations were performed. The 
first was carried out under equilibrium NVE conditions and the force on one of the 
nanotubes (Figure 4) was calculated in the direction normal to the nanotube 
surface. In the second simulation the same nanotube was loaded incrementally with 
axial force on each atom until the nanotube began to slide. Both the velocity of the 
nanotube and the distance that it traveled were monitored throughout the 
simulation. Two definitions were used to determine the initiation of nanotube axial 
motion. In the first the nanotube velocity was consistently positive above the 
thermal induced velocity fluctuations. In the second definition the nanotube was 
considered to be moving when it had traveled for 0.2 nm. The friction coefficient is 
calculated as the negative of the ratio of the applied to the normal load and is 
presented in Table 2. As explained in the following discussion. The normal load 
and apparent friction between nanotubes in a bundle may prove to serve as a 
starting point for macro-scale calculations of yam fiber friction, a quantity that 
plays a role in both strength and stiffness calculations. 



Table 1 . Stress on the Nanotube System under Compression 


Percent 

Compression 

(Nominal) 

Stress Component in 
Transverse or x- 
direction (GPa) 

Stress Component in 
Axial or z-direction 
(GPa) 

0.0 

-0.18 

-0.04 

2.0 

1.73 

0.80 

5.0 

1.77 

0.97 

10.0 

1.86 

1.24 

15.0 

1.89 

1.64 


Table 2. Apparent Friction Coefficient for Nanotubes Sliding. 


Percent 

Compression 

Method of 
Determining 
Pull-Through 
Force 

Pull-Through 

Force 

(pN/atom) 

Normal 

Force 

(pN/atom) 

Apparent 

Friction 

Coefficient 






0.0 

velocity 

0.053 

-0.277 

0.190 


distance 

0.065 


0.230 






2.0 

velocity 

0.120 

-8.40 

0.015 


distance 

0.120 


0.014 






5.0 

velocity 

0.075 

-9.00 

0.0083 


distance 

0.068 


0.0075 






10.0 

velocity 

0.083 

-5.40 

0.015 


distance 

0.081 


0.015 






15.0 

velocity 

0.043 

-4.930 

0.0088 


distance 

0.030 


0.0060 


During compression of the nanotubes system, the relationship between the 
stress in the transverse direction and the normal load on the nanotube is non-linear. 
This non-linearity is not surprising because the nanotubes deform under 
compression. 






Figure 4. Molecular structures of carbon nanotubes simulated with and with out transverse 
compression, (a) zero compression (b) 5% and (c) 15% compression. The marked nanotube 
(x) was loaded in the axial direction. 

The calculated friction coefficient decreases relative to the undeformed 
system as the system is compressed. The primary reason for this is that the normal 
force on the nanotube is also decreasing. The pull-out force initially increases at 
2% compression then reduces with increasing compression. It is expected that 
during scale-up to the mechanics methods that the normal force on the nanotubes 
could be related to lateral forces induced under twist and axial load of a single yam. 

The polymer-yam-nanotube interfaces have yet to be fully explored so it is 
worth noting that a nanotube sliding on a non-smooth surface such as the inside of 
a functionalized outer nanotube shows different friction characteristics than when it 
slides on a smooth surface. When the outer nanotube is not functionalized, the 
inner one slips easily in the axial direction. Therefore, the nanotube has a low 
friction coefficient in the non-functionalized case. When the outer nanotube is 
functionalized, the (absolute value) magnitude of the normal force is much smaller 
than the non-functionalized case, and inner one ‘sticks’ while trying to move 
axially within the rough surface of the outer nanotube. Both of these changes in the 
nanotube physical behavior raise the nanotube friction by at least two orders of 
magnitude. As demonstrated by others, [14], additional changes to nanotube sliding 
behavior occurs when there are defects in the outer nanotube or covalent bonding 
between nanotubes in MWCNTs. As described by Frankland [13], investigations of 
the composite should also consider covalent bonding of the nanotubes to the 
polymer. 

At the atomistic scale, MD simulations by Qian [15] were used to estimate 
the load transferred to the central nanotube as a nanotube bundle of seven 
nanotubes is twisted. The simulation calculates the transferred load in the bundle 
versus the twist angle from 0-210 degrees. At zero degrees, a force of 0.048 eV/A 
is transferred. At 120 degree the transferred load peaks at 1.63 eV/A. For the 
current MWCNT yam structure, the degree of twist of nanotubes at the nanoscale is 
expected to be small relative to the overall yam twist. For the twisted bundle of 
fibers, it is assumed that the behavior of interest is directly related to the yam twist 
and load-induced lateral compression that occurs at several length scales. 



Macro-level analysis 

At the continuum mechanics scale, a large part of the relevant textile 
literature deals with dry single yams, i.e. a yarn without any supporting matrix. A 
short review of several papers that describe dry yam stress-strain behavior follows. 

The complex geometry associated with twist and material variations implies 
that modulus of a yam is a volume averaged structural property rather than a 
material property. This complexity results in both structure and axial load levels 
affecting all the mechanical properties of a twisted yam. Generally, yam twist is 
characterized by a twist multiplier (TM) [9]. 


TM = C^ 2 t 


( 1 ) 


Where r is the number of twists per unit length and C is the linear density 
with units of mass per unit length. During tensile loading of the yarn, the apparent 
stress-strain behavior will generally have a short linear region characterized by 
slippage of the fibers followed by nonlinearity with considerable extension while 
the fibers continue to slip and finally fiber failure occurs. During this axial loading 
of the yam, off-angle forces in the yam are generated by yam twist. This multiaxial 
load state consists of combinations of axial tension and lateral compression. 

Rao and Farris [16] considered the effects of material anisotropy on twisted 
yams. They assumed that the yam was composed of a series of thin-walled layers 
that resembles a unidirectional composite. The layers are located at a radial position 
(r) in a yam of radius ( R ). The location of the filament is further defined by the 
angle (0) between the single filament at r and the yarn major axis. Further, angle 
(a) is the angle (0) at radius R. The apparent Young’s modulus of a twisted yam is 
calculated from the following expression 

E = — - — f — tan# sec 2 6d0 (2) 

tan 2 a } ° S n 


where Sn is the principal compliance in the plane of a single transversely isotropic 
layer. If one assumes low twist and a < 15° then Rao arrives at an expression for 
the apparent Young’s modulus of the yam. 


E = E. 


A 3r o +P 

v 2dT o J 
d = 

E. 


, (1-0 2 ln (i -d)r 0+ d 

<f 3 tan 2 a T 


T = cos a 


( 3 ) 


The terms E z and E s are the longitudinal and shear modulus respectively of the 
fibers. 




Figure 5. Schematic representation of relationship of coherence, obliquity, and twist to yarn 
strength. 

Yam strength has been found to be a function coherence (the degree of fiber 
packing) and obliquity (degree of order or parallelism of the fibers in the yam. 
These basic relationships are visualized in the following schematic representation. 

Considering yam strength, there is an “optimum” twist associated with 
maximum strength of the yam. Below this optimum, the yam will fail due to fiber 
slippage while above the optimum; the yam will fail by fiber fracture. As presented 
by Pan [17], [18] as a consequence of the lateral compression forces on the yam, 
the yarn behaves as a chain of independent fragments. Each fragment is an 
individual due to fiber non-uniformity and specific fragment length. During yam 
extension, fiber breakage continues until the length of the fiber fragment reaches a 
minimum dictated by the load carrying capability, or strength, of the fragment. This 
minimum length can be considered the “critical length” for this structure. The 
critical length has also been shown to decrease with increasing lateral pressure. 
Twist of the fiber and yam diameter also plays a role in determining critical length 
although it is thought that twist improves the strength of high-performance fibers 
but decreases the strength of low-performance fibers [16]. 

Based on statistical interpretation of experimental data, Pan has proposed 
the following strength relationship. 



f i \ 
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( 4 ) 


Where cr y is the yam strength, <y h is the fiber bundle strength, l f , / . are the 
fiber length and critical length respectively, v f is the volume fraction of packing 
density of the yam, and /?,7 are the Weibull shape parameter and orientation 



efficiency factor respectively. An expression for the critical length of a yam was 
presented by Realff [19] 


2T, 

i = L 

71 d JLlCJ l 


( 5 ) 


Where, 7/ is the fiber tensile force at yam fracture, d is the yam diameter, cti 
is the lateral pressure, and p is a frictional coefficient. Realff [19] also presented 
this expression but asserted that contrary to Hearle [20], twist induced lateral 
pressure raises the overall strength of a yam by increasing the apparent strength of 
the fibers. 

Investigations by Gosh [21], [22] have demonstrated yam mechanical 
behavior for a variety of material systems. Their studies have highlighted the 
influence of the helical geometry of the twisted yarn and the role of the resultant 
multiaxial forces generated during yam axial tension by the unique structural 
arrangement of the fibers. Based upon test observations, they have proposed the 
following empirical expression for yam tenacity. 


0 = 



— cos 2 0 
100 


( 6 ) 


Where the terms of equation (6), defined below, are determined from test 
and observation. The concepts of yam tenacity and the relationships to stress can be 
found in [20]. 

n h = number of fibers at location of yam break 
n e = average number of fibers in yam cross-section 
n h = n e -1.96(7 

<y = standard deviation for 95% confidence level of number of fibers in yam cross-section 
F h = fiber bundle tenacity at gage length 
6 = average helix angle of fibers in yam failure zone 
3 b = percent of broken fibers in yarn failure zone 


Experimental evaluation of the nanotubes-polymer composite, described in 
the Materials section, reveals yam stress-strain behavior in the composite similar to 
stress-strain behavior of a single dry fiber. The four-yam polymer composite used 
in this study was subjected to uniaxial tension in a test stand with load measured 
using a standard load cell calibrated to +/-10N. Stress is engineering stress based on 
the load cell data and the average specimen cross-sectional area. The static tests are 
conducted with a small tensile preload and a constant displacement rate of 
17mm/min for all the specimens. Once loaded to 50% of the gage length, the 
sample is unloaded at the same rate. Axial strain is calculated as the change in grip 
or crosshead displacement divided by the initial crosshead displacement as 
measured by the test machine’s displacement transducer. Load and displacement 
data are collected at a rate of 2 samples/sec during the test using test machine’s 



digital data acquisition and control system. After the static tests were completed, 
the specimens are examined for local yarn failure characteristics such as location, 
mode and yam-to-polymer interactions. Strain to failure levels in the polymer 
exceeded the designed limits of the test and consequently none of the composite 
specimens failed globally. In order to identify features at the nanoscale (type of 
fracture, interface characteristics, etc.) that will help provide insight into the failure 
of the yam composite, optical and electron microscopy is performed after testing. 
Figure 6 provides some sample images of a cut and fractured yam-polymer surface. 
Note the presence of polymer matrix in-between yams in image 1 and the pull-out 
of the fractured yam away from the surrounding polymer in images 2 and 3. 



Figure 6. SEM images of a four ply yarn. Image 1: cut surface [23]. Image 2: Fractured 
surface. Image 3: Enlargement of fractured yarn. 

A representative stress-strain curve for the static tensile test of the 
composite is provided in Figure 7. Upon loading, initial failure corresponds to the 
maximum stress (point 3). At initial failure, all five yams failed simultaneously at 
the same location with respect to the longitudinal axis (Figure 8). The majority of 
test specimens had initial failure occur in the gage section. Upon further loading, 
the specimen has lost stiffness but continues to carry load until the stress reaches 
the secondary failure level (point 5). At this point simultaneous failure of all five 
yams occurs again but at a location different from the point of initial failure. As the 
yam fails in distinct segments, the effective stiffness of the composite trends 
toward the stiffness of the polymer matrix. Finally, upon unloading, the specimen 
does not return to its initial configuration (point 8). 
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Figure 7. Representative Engineering Stress - Strain curve for 5 yarn composite loaded in 
static tension. Note: 1 - 1% strain; 2 - 4.5% strain; 3 - first failure across all yarn; 4 - 
8%strain; 5 - secondary failure; 6 - 20%strain; 7 - 50% gage length (point of unload); 8 - 
final strain 

At failure, the yarns did not pull-out from the matrix but did exhibit spring- 
back or recoil along the loading axis. This behavior is highlighted in Figure 8 by a 
series of photographs of a composite after failure. The failed yams take on a 
resemblance to a buckled column structure on an elastic foundation. This type of 
behavior was also observed in [6], Their explanation revolves around the model 
that predicts shear and compressive loads along the interface between yarn and 
matrix. These loads are due to differences in yam and matrix longitudinal strains 
and lateral contractions. Yam failure in Figure 7 was clearly defined and the post- 
mortem inspections revealed that the yam breaks were even with each yam single 
failing on the same cross-sectional plane. As illustrated in Figure 9, the buckled 
fiber contributes little to additional load carrying capacity of the composite. In 
addition, examining Figure 9, it is reasonable to contribute the regularity of the 
breaks and yam segment lengths to the ideas behind the concepts of critical length 
as presented in equation(5). 





Figure 8. Optical micrograph of failed tensile specimen gage section-longitudinal view. 12x - 
view of failure zone of all 5 yarns; 50x - inset of yarn 1 and yarn 2; lOOx - inset of yarn 1; 
150x - inset of yarn 1 and initial configuration of yarn 1. Note: Vertical red line denotes the 
center of the sample. 
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Figure 9. Optical microphotograph (12x) of failed composite tensile specimen illustrating the 
sequence of yarn segments. 




SUMMARY 


A novel nanostructured composite material was considered as a candidate 
for multi-scale analysis of strength and stiffness. As discussed by [24] and [25], it 
is often the interface region that governs the material response. This is true for a 
range of properties including mechanical, electrical and thermal. The implication is 
that the multiscale analysis method must provide sufficient detail in the interface. 
Due to practical limitations on computational speed and size, the level of detail 
must be tempered against expediency and efficiency. The material interface 
requires a transition between scales and their associated methods. This paper 
presented an initial attempt to experimentally and analytically describe stiffness and 
strength of a MWCNT yam-polymer composite. 

The matrix material in the composite is an elastomer and the reinforcement 
is provided by five, four-ply multi-wall carbon nanotube yams. The composite has 
a yam volume fraction of 0.0028. An experimental study was performed to 
determine the static mechanical response of the multiwall carbon nanotube - 
polymer composite. Strength, stiffness, and failure were characterized for both the 
polymer and the composite loaded in uniaxial tension. The yam response is 
dependent upon inherent material properties as well as yam structure, geometry and 
manufacturing method. Despite the low yam volume fraction, the composite 
modulus was increased by over a factor of ten compared to the base polymer. The 
composite exhibited linear elastic behavior up to initial failure. The high strain 
value (~5%) at failure resulted in catastrophic yarn failure with all yams failing 
simultaneously at some characteristic length. Failure modes of the yams were not 
unlike fabric yam failure modes described in the textile literature. As is typical with 
twisted yams, the final yam structure or bundle does not retain the strength or 
stiffness of the base material or strands. The scale-up from MWCNT to the yam 
implies the final yam properties are a function of base material properties, material 
twisting, fiber migration, frictional forces and lateral compressive forces developed 
at each material level or length scale. 

Due to the homogenization process, a set of assumptions accompanies each 
scale in the yam structure including the interface transition regions. The solutions 
to multiscale analysis are material and application specific and transition between 
scales is the essential challenge. For a novel material such as the multi-wall carbon 
nanotube yam, the solution to predicting mechanical behavior requires the 
coordinated contributions of engineering, science and computational expertise as 
well as the direct coupling of experimentation and modeling. 
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